CO 

CD 

o 

C3 

<c 


STABILITY  OF  THE  FOUR  IER  METHOD 


by 


Heinz-Otto  Kreiss  and  Joseph  Oliger 


STAN-CS-77-616 
AUGUST  1977 


COMPUTER  SCIENCE  DEPARTMENT 
School  of  Humanities  and  Sciences 
STANFORD  UNIVERSITY 


D D C 


L 


Unclassified  

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whan  Data  I'.na.-a J) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1 . REPOS T NUMBER 

STAN-CS -77-016 

2 GOVT  ACCESSION  NO. 

3.  REClPl-fcNlLii  CATALOG  NUMBER 

UT.’-'  ' • - . 

4.  TITLE  (and  Subtitle)  

5.  TYPE  OF  REPORT  4 PERIOD  CdtEREO 

STABILITY  OF  THE  FOURIER  METHOD, 

-r  t *r 

Technical,  August  1977 

JVs' 

OPG.  REPORT  NUMBER 

SIAH-CS-77-olo  A 

7.  AU  THORfs) 

' Heinz -Otto I Kreiss  Mfie l Joseph/ oiiger  1 

8.  CONTRACT  OR  GRANT  NUMBER^ 

on  if  Nd!doi4-7b  -C-1132J' 

* — 1 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Stanford  University 
'Computer  Science  Department 

Stanford,  Ca.  oL30b 

10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  A WORK  UNIT  NUMBERS 

, 

/ > y ' + , 

11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Office  of  Naval  Research 

j'{\ 

12.  REPORT  DATE 

Department  of  the  Navy 

Arlington,  .a.  22217 

13.  NUMBER  OF  PAGES 

28 

l«  MONITORING  AGENCY  NAME  4 ADDRESS!-!!  dt  Iterant  Iron.  Controlling  Office) 

01'..-  Representative:  Phillip  S’irra 

Durand  Aeronautics  Bldr.,  : . lob 

15.  SECURITY  CLASS,  (ol  thin  report) 

Unclassified 

Stanford  Iniversity 

1 Staufor  i,  Ja.  'i4?0b 

1 5a.  DEC  L ASSIFI  CATION/ DOWNGRADING 
SCHEDULE 

15  DISTRIBUTION  STATEMENT  (of  this  Report) 


Keloasable  without  limitations  on  dissemination 


'distribution  statO/ii:t  a 

Approved  for  public  release; 

Distribution  Unlimited 


17  DISTRIBU^,  _'N  STATEMENT  (of  the  abstract  entered  in  Block  20.  if  different  from  Report) 


13  S'JPPL  EMEnt  ARY  NOTES 


19  KF.Y  *ORDS  (Continue  on  reverse  side  if  necessary  and  identify  by  block  number) 


20  ADSTRACT  ( Continue  on  reverse  .ere/-  if  necessary  and  Identify  by  block  number) 

— In  this  paper  develop' a stability  theory  for  the  Fourier  (or  pseudo- 
spectral)  method  for  linear  hyperbolic  and  parabolic  partial  differential 
equations  with  variable  coefficients. 


Unclassified 


DD  1 j AN  ^73  1473  EDITION  OF  I NOV  65  |5  OBSOLETE 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (W>«n  D ala  Enlara^) 


STABILITY  OF  THE  "OURIER  METHOD 


* 

Heinz-Otto  Kreiss 


and  Joseph  Oliger 


ABSTRACT 

In  this  paper  we  develop  a stability  theory  for  the  Fourier 
vor  pseudo-spectral  method  for  linear  hyperbolic  and  parabolic  partial 
differential  equations  with  variable  coefficients. 


Courant  instttute  of  Mathematical  Sciences,  Hew  ork  University, 
251  Mercer  Street,  N.Y.  10012.  Supported  in  part  by  the  National 
Science  Foundation  under  Grant  ATM  76-10218. 

** 

Department  of  Computer  Science,  Stanford  University,  Stanford,  CA 
1A305.  Supported  in  part  by  the  Office  of  Naval  Research  under 
Contract  N0001L-75-C-1132. 


DISTRIBUTION  STATEMENT  A 

Approved  tor  public  release; 
Distribution  Unlimited 


A 


1. 


Introduction 


The  collocation  method  based  on  trigonometric  interpolation  is 
called  the  Fourier  (or  pseudo- spactral)  method.  It  has  been  used 
extensively  for  the  computation  of  approximate  solutions  of  partial 
differential  equations  with  periodic  solutions.  A satisfactory 
theoretical  justification  for  equations  with  variable  coefficients 
has  only  existed  for  equations  written  in  skew  symmetric  form  [3 , 6,  ?]• 
Recent  work  of  Majda,  McDonough  and  Osher  [8]  treats  hyperbolic  systems 
with  C°°  coefficients. 

In  this  paper  we  develop  a stability  theory  for  linear  hyperbolic 
and  parabolic  partial  differential  equations  with  variable  coefficients. 
The  generalization  of  these  results  to  nonlinear  equations  follows  if  the 
problem  has  a sufficiently  smooth  solution.  We  restrict  our  discussion 
to  problems  in  one  space  dimension.  The  extension  to  problems  in  more 
space  dimensions  is  immediate.  Error  estimates  can  easily  be  derived 
using  our  results  following  those  in  Kreiss  and  Oliger  [7]  and  Fomberg 

[5]  • 


. ..  .'on  □ 
C 


CISOT-.Y  ;•  ;$ 


2.  Trigonometric  Interpolation 

In  this  section  we  collect  s one  known  results  on  trigonometric 
interpolation  (see  [ k,  5,  7^-  Let  N be  a natural  number, 
h - (2N+1)”1,  and.  define  grid  points  xy  - vh,  v - 0, 1,2, . . . ,2N. 
Consider  a one-periodic  function  v(x  , vix.  •-  v(xfl),  whose  values 


vv  - v(x  ) are  kn°wn  at  the  gridpoints  Xy.  We  define  a discrete 
scalar  product  and  norm  by 


2N 


(2.1)  Uux),v(x))  = I u-.x  v(x  "'h  , 

v=o 


(u,u  h 


The  trigonometric  polynomial  w(x 1 of  degree  I!  which  interpolate 


v ,x ) in  the  points  xy,  i.e. 


(2.2! 


w(x,,)  = v(x. 


V . , 1,2 , . . . ,2ii  \ 


is  uniquely  given  by 


(2. J) 


wvx)  = a (o  )e 

ai=--N 


2iriux 


where 


2.U) 


i \ \ 2TTicux . 

a(<n)  = (v(x/,e  ), 


This  follows  from  the  orthonormality  of  the  exponential  function, 

/ 

0 if  0 < | m-n | < 2N 


(2.3) 


( 


27Tinx  2TTimx 


V- 


1 if  m n . 


The  usefulness  of  trigonometric  interpolation  stems  from  the  fact 
that  the  smoothness  properties  of  the  function  are  preserved  and  that 
the  convergence  is  rapid  for  smooth  functions.  Let  the  L0 -scalar 

d. 

product  and  norm  be  defined  by 

(2.6)  (u,v)  = /*  uvdx  , ||u|j2  = (u,u) 

0 

We  will  need  the  following  well  known  theorem. 

Theorem  2.1.  If  w^,w2  interpolate  v^  and  v^,  respectively,  then 

(2.7)  (Wp’Wp )h  = Vw1,w2 ) = (v  ,v  ) and 

12.8)  !|wl(x)!|2  = ||vL(x)]^  = Z | a (00 ) |2 

u)=-N 

It  will  be  convenient  to  work  with  the  following  class  of  functions. 

Definition  2.1.  P(a>M)  is  the  class  of  all  functions  v(x)  which 
can  be  developed  in  a Fourier  series 

(2.9)  v(x)  = Z v(cu)e27rimx 

03=  -00 

with 


(2.10)  7'.  I [ |27roo|a  * l]v(a))  I2  < m2  . 

OD--00 

P(a>M)  is  contained  in  the  Sobelev  space  H^. 

We  now  need  the  relationship  between  the  Fourier  coefficients  v(a>) 


k 


of  a given  function  vix)  and  the  coefficients  al“)  of  its 
trigonometric  interpolant  wlx).  This  is  contained  in  the  following 
well  known  result  [ 4,7]. 

heorem  Let  v be  given  by  12.9'  and  w given  by  (2-5 i and 

(2 .4  1 then 


v2.ll'  aiu>)  = ’ v (cn-j  1211  *-1 1 ' , I co | < j\j 

j:  -OO 

We  can  now  investigate  the  rale  of  convergence  of  the  interpolating 
polynomial  to  a function  vlx  Pla,M). 

oreri  ■ • . Let  wx'  Pva,M)  with  y > l/2.  Then 


.12) 


V X -W  \X  ' < 


127 m 


127711  ) 


2a 


z — — 

j=-l  12,7 -l)1- 


1 1/2 


a 


127m )' 


a 


1 

where  ' 1 *■  2 

a j=l  12J-1)2« 


roof.  We  write  (2.9)  as  vix'  v,. ix)  + v.  ix)  where 
“ r.  R 


N ^ . 

v ri  ~ i 27no3x  x v - \ 27ricox 

vTIvx»  2:  v(a>)e  , v ix)  = L vVa)e 

a--N  M>N 


<e‘  w„.ix’  and  w.,  (x 1 be  the  trigonometric  interpolants  of  v,j(x) 
and  v.  x , respectively.  They  are  given  by 


N 


w.,  lx)  = £ a'  loo  ;e 


IN)  ^)_27ria>x  _ (N ) 


lav)  - (v.Iix'l,etTicuXV 


\ v (R ) \ 2iriu)x  (R ) / /„  (v\  2Triujx^ 

w lx)  = Z a (U5)e  , a (u> ) = (,vR  vx ) , e 


The  trigonometric  interpolant  of  v(x)  is 


w(x)  = WR(x)  + wR(x) 


w (x)  interpolates  vR(x)  in  the  2N+1  points  of  (2.2),  and  from 
{2.3)  we  have 

wN(x)  = vN(x)  . 


Therefore, 


jv(x)-w(,x)!!2  = ||vr(x)-wr(xH|2  = ltvRlx)l!2  + )lwR(x)l|c 


since  vD(x)  is  orthogonal  to  wR(x'.  By  (2.10)  we  can  write 


- (CJD ) = 


| pmo  | a + 1 


where 


Z |vla>)|2<,vp 


Therefore , 


!fvo(x)|P  - z |v(co)  |: 
R hi  > N 


> N |2T7tu|a+l 


^ 


By  Theorem  2.2 


6 


u 


N ,nv 
v,  ta'"' 
cu- -N 


N 00 

I I Z vi03+j  (2N+1>  I !‘ 

03-  -K  j=-co 

it  0 


N 

£ 

03= -N 


CO 

vyo3+.)  (2H*-1 ' 1 

j---®  |2TT(0'-t-j  vcN*-l)  ' |°+1 

j/o 


N 

< z 


03= -M 


— ~ p*  v-  I v(ov»-j  (2N+1 ' ) )2 

j = -00  { |27T(0>t-j  ( M + l " |Q+  \Y  «j--co 

jto  j/o 


2Rp 


jrt\f 


\2  a 


Z (23 -l) 

j=i 


a 


a:..  -orem  follows. 


Remark.  >bserve  that  the  contributions  to  the  error  by  v and  w 
• ri  R 

are  of  the  same  order  if  a , 1/2.  v is  often  called  the  aliasing 
error.  Thus,  we  see  ‘hat  if  v is  at  all  smooth,  then  aliasing 
rlays  no  important  role. 


"he  following  result  follows  immediately  from  the  last  theorem. 

’urollary  " ■ 1.  Let  v(x ' Plco'4'  with  a ,i  *-  1/2  , ,)  a natural 
number.  Then 


2.1/ 


4v<*> 

dx'5 


MC 


(2im)a~5 


3.  Stability  of  Fourier  Methods 

Let  vlx^  be  a one-periodic  function  whose  values.  v(xy  are 
known  at  the  gridpoints  xy  = vh,  h = (2N+l)  1.  If  we  want  to  approxi- 
mate dv(xv)/dx  we  can  compute  the  trigonometric  interpolant  (2.3'1  of 
v(x',  differentiate  it,  and  use  its  derivative 

N 2TTicuxv 

(3.1)  dw (.x  )/dx  = Z (2TTia')a(o>)e 

V ut-N 

as  an  approximation  of  dv(xy)/dx.  The  computation  of  (.3-1'  in  all 
of  the  gridpoints  x , v = 0,1,2, ... ,2N  can  be  done  using  two  discrete 
Fourier  transforms  vDFT^  and  2N  complex  multiplications.  Also,  if 
we  know  that  v(x)  with  a>  3/2,  then  Corollary  2.1  gives 

us  the  error  estimate 

MC  , 

(3.2)  dv/dx  - dw/dx  < - — j . 

(.2Tffl)n 

Higher  derivatives  can  be  computed  analogously. 

The  above  process  is  linear  so  it  can  also  be  represented  using 
matrix  notation.  Let 

v = (v(xQ),. . .jvlxgjj))7  , y = vdw  (,xQ ) /dx , . . . , dw  kx^ ) /dx ) ' 1 

denote  the  (2N+1)  dimensional  vector  formed  of  the  grid  values  of 
v{x)  and  dw/dx,  respectively.  Then  there  is  a 12N+1)  X (2N+1) 

If  y is  a vector  then  y denotes  its  transpose  and  y its 
conjugate  transpose.  The  same  notation  will  be  used  for  matrices. 
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where 


Ou  - i-vbdu/dx  + (Ubu)/dx)  , Ru  = - | db/dx  u . 

Q is  skew  Hermitian  and  R is  bounded.  There  are  many  problems  where 
R = 0.  For  example,  we  can  write  udu/dx  in  the  form 

1 2 1 
udu/dx  = --(udu/dx  + du^/dxl 

!«cw  consider  the  partial  differential  equation 


u - bvx)bu/bx  = Qu  *-  Ru  , u = bu/bt  , 


then 


u,u)  - ,u,u, ' + vu, ,u^  - vU,Qu'  + lQu,u)  + vu,Ru'  + vKu,u<  -(u,udb/dx 

t t. 

ar.d  we  have  an  energy  estimate.  If  we  approximate  the  above  problem  by 


dv 


1 i l 

— - -vbS  + Sb)v  - ^ b v 


dt  2 


where 


b lx  ' 
o 


b (x  ' C ...  0 


0 b lx, 


2N 


1 x- 

db  ^x  ^ 
o 

n . . . 

0 

dx 

0 

db  (x^ ) 

0 . . . o 

dx 

db  ^x. 

2N 


dx 


Then  we  obtain  the  same  energy  estimate  because 

(bS  + SfJ)  = - ISS  + Sft) 


10 


is  skew  'lermiti an  and  therefore 


- (v,F  v) 
— 7 x— 


h 


F!ie  above  procedure  can  be  generalized  considerably.  Consider  the 
parabolic  system 


l,3.?’1  u = iAu  ) + Bu  +-  Cu  , u = 6u/bx 

t XX  X X 

where  u denotes  a vector  function  with  n components,  A,  B,  and  C 

are  n X n matrices,  A and  B are  Hermitian,  A is  positive  definite,  and 

C and  dB/dx  are  uniformly  bounded.  We  can  rewrite  this  system  in  the  form 

(?.6)  u,  = (,Au  ) + — (B  u +■  (B  u)  ) + C u , 

t x x 2 x x 1 

where 

C1  = C - i bB/dx  . 

We  then  obtain  the  energy  estimate 

(u,u).  = -2(u  ,A  u ) + 2 Real  vu,C.u' 

t X X 1 

which  depends  solely  on  the  property  that  6/5x  is  skew  Hermitian. 

Thus,  we  obtain  a corresponding  estimate  if  we  replace  6/dx  by  S 
and  approximate  (J.6)  by 

dv , _ „ 

'I/O  — = 5 ASv  » S + SB  )v  + C v 

Q.L>  c.  X 


li 


The  estimate  is 


d_ 

dt 


+ C,  )v). 


J1  '-'h 


where  we  extend  our  earlier  definitions  of  the  discrete  norm  and  inner  pro- 
duct in  the  obvious  way.  Here  v is  the  vector  with  vector  coir^onents 
v(x^)  and  A,  B,  C.. , and  S are  block  diagonal  matrices  with  blocks 
A(x^),  B(x^),  C^(x,  ),  and  S,  respectively. 

The  system  of  ordinary  differential  equations  (>3*7)  can  be  solved 
using  an  appropriate  difference  method  for  ordinary  differential 
equations.  However,  the  approximation  (3-7)  requires  about  twice  as 
much  work  as  the  simpler  approximation 


(3-8) 


dv 

— = SASv  + BSv  + Cv 
dt  — — — 


of  (3.5).  Since  numerical  experience  has  shown  that  approximations 
of  the  form  (3.8)  can  be  unstable,  it  is  desirable  to  find  ways  of 
stabilizing  them  which  are  cheaper  to  use  than  reverting  to  (3.7).  We 
can  achieve  this  by  adding  appropriate  dissipative  or  projective 
operators.  We  will  now  develop  this  approach  in  deta’l. 

It  is  easier  to  do  this  if  we  work  within  the  space  T of 

INI 

trigonometric  polynomials 

(3*9)  p(x)  = i.  p(u>)e277^ax 

A vector  function  v(x)  or  a matrix  function  B(x^  will  belong  to 
Tjj  if  all  their  components  do.  There  is  a one-to-one  correspondence 
between  a polynomial  13*9)  and  its  values 
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V - Mxq),...,v(x2N))'  . 

Thus,  there  is  a linear  operator  P such  that 

Pv  = v(x)  , i.e.,  Vlxv^  - Vy  , V = 0,1,2,. 

If  vlx"1  Tjj  then 


v 3 . 10 ) PSv  = dv/dx  . 


Let  Six'1.  v(x^  t T...  Then  we  define  w(,x)  = B(x)*v(x) 

IN 

convolution 


0.1l)  wi,x^  = BvxVv\x'>  = Z w(v)e2TT^Vx 

V=-N 

with 

/ 

+N 

£ B(u)iv\.v-a)  +■  v(V-2N-l~n))  for 

u=-N 

(3-12>  wvV)  = < 

N 

£ B(u)(v(v-al  + v(v+2N+l-ii ) ) for 


where  we  have  used  the  convention  that  v(<*>)  = B(a>)  = 0 
BvxVix'  is  a trigonometric  polynomial  of  order  2N. 
its  interpolant  is  given  by  B(x)*v(x).  Therefore, 


3-131  w(x)  = P(Bv)  - Blx)*v(x) 


. . ,2N  . 


to  be  the 


V > 0 

V < 0 . 

if  M > N. 
By  theorem  2.2 
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Lemma  3*2.  Let  B(x)  e T„,  be  a matrix  and  v.w  e T„  be  vector  functions. 
N N 

Then 

|(v,B*v)|  < max  |B(x)|  • ||w|(  (|v||  . 

“ ° < x < 1 

and,  if  B is  Hermit ian, 

(w,B*v)  = (B*w,v)  . 

Proof.  By  theorem  2.1  and  (3*13) 

(w,B*v)  = (w,B*v)^  = ^WjBv)^ 

If  B is  Hermit ian,  then 

(w,iv)h  = (Bw,v)h  = (B*-u,v)h  = (B*v,v)  . 

Also, 

| (v, Bv ) I < |B|||w|| h||v||  = max  |B(x  )|||v||||v||  • 

n 0 < xt  < 1 

and  the  lemma  is  proved. 

We  can  now  wr^te  equation  0*8)  as  an  evolution  equation  in 
via  the  isomorphiam  P. 

(3-lM  vt  = (Vvx}x  f Vvx  * CN*V 

where  ,B^  ,Cj.  and  v are  the  trigonometric  polynomials  in  T,, 
which  interpolate  the  discrete  values  A(xv),  Blxv),  Ctx 

respectively.  The  term  w = B.*v  can  be  written  as 

N X 

w = Bn*vx  = Qy  + Rv 

1U 


whei'e 


13-15) 


**  ■ 5'b»*vx  • 'VvV 


Kv  - -U<  kv  - B *v'  ) 

1.  x N x 


It  follows  from  lemma  3-2  that  the  operator  q is  skew  Hermitian. 
straight  forward  application  of  ■,  . 11  gives  us 


Rv  = R1  v + R v , 

-L 


?.  iu)e2iriu* 


where 


(3.1o)  r^(co)  = _rri 


^3-17)  r,  (u>)  _ yri  (2N *-1 ' 


pBN(u)vv^-.l  * v^«'-?N-l-.. ) ' for  u>  ; 0 


uBnU  AV  ~l!  ) + v(um-2N+1-.,  ' for  O)  < 0 


BfI  )v(o.'-2N-l-u ) for  cu  p 0 


Bp. (a ) v(a)+2N+l-|_i ) for  a)  < ( 

n--N 


By  13.12) 


(3.16) 


R1V  " 2 dVdx  * 


Therefore,  by  lemma  3.2,  the  operator  is  bounded  if  B £ P(q,M)  with 

' > ->l?>  certainly  if  B is  twice  continuously  differentiable  (see  [1]). 

In  general  we  can  not  expect  that  (v,R^v)  is  bounded  independent 
of  N.  For  example,  if  flue)  - 1(1  ► i sin  2irx)  then 
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b^o)  = i,  sNa)  - -V_1)  = - t lf  ®n(u))  = 0 if  1“!  + °’1 


and 


f2(N)  = £ (2N+l)v(-N)  , r2(-N)  = £ (2N+l)v(N)  , r^)  =0  if  M f N . 


Therefore,  by  Parseval's  relation. 


iv,R,v)  = | (2N+1)  Real{v(N)v(-N)}  . 


Now  assume  that  there  are  constants  and  fi  > 1,  independent  of  N,  such 
that 

M 

(M9)  IB  U ) | < ~r  for  n t 0 . 

|27tJP 


Then  we  obtain 


N N 

| \v,R0vH  < tt(2N*-1)(|  E E BN(u)v(a>-2N-l-u)  I 

" “ O3=o  u=-N 

*.2C)  

-1  N 

- I E v (o3 ) E B,Eu  )v(om-2N+1-u  ) | ) 

ufc-N  u— N 1 


wV '>re  v - - 0 for  |t|  > N.  By  (3-19^ 


1 v \U) ' E B.  U)v(,co-2N-l-n)|  < 

*0  U--N  r’ 


— - E |v(o>)|  |v(^-2N-1-u)|  < 

1 „ -N  \*,J»  «*0 
M/*0 

M.  E • E | v (o) ) | |v(o>-2N-l-p.)  | < 

a=-N  |27Ta|P  oteN+u+l 
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<«•- 


M,  -1  , 

1 v 1 


-N 


|27TU|P 


N 

v I ~ I -J  I " . , 

liviu  I *•  ' V -1- 


M1  -V1  1 

L, 


“ u=-K  |27Tu!P 


v'vv:'|'  f !v(,-u>)|‘M 


< 


N ^ N 

7T  - lv^)|2  Z 


ufc-N 

<V0 


i=N-|o)|+-l  ^2-a|P 


rhere  is  a constant  such  that 


z 


*1 


i=N-  | so | +1  (?7T.. 'P  (N-  |cd|+1  ) 


6-1  ' 


Ki  (V2rr)  ( 5/(  r -1 ) ) will  do.  Furthermore,  the  same  estimate  holds  for 
the  second  sum  on  the  right  side  of  (3.20).  We  obtain 


3.;  i 


where 


i>v. 


R2v  1 ! < M1K1  • Z rjvi^) |2 


o^-N 


i2N+l'TT 

ra>  = — — — — ^ “ f 0 , T = 0 

(N-|a)  Ul)“  1 ° 


Consider  the  system  (3.1k).  We  have,  using  (J.15)  and  (3.16), 


i v, v ) , = 2 Real  { (v,  (A^v^ ) t-  lv,Qv)  + (v,Rv)  + (v,C  *v)} 


J.22) 


‘^VV^x’  f 2 Real  (v^cn  ‘ | &BN/&x)  * v)  +•  2 Real  ^v,R0v 
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A is  positive  definite  by  assumption,  i.e.,  there  is  a constant 


a > 0 such  that  A > cl*  Therefore, 


(v  ,A..*v  ) = (v  ,Av  ),  > a'!v 
x 71  x x x h - 1 x 


By  Parseval's  relation  and  (3»2l) 


-2  (v  ,A_r*v  ) +•  2 Real  (v,R_v)  < 
x 71  x 2 — 


0.2?: 


2 S (-o-(2nco)2+  M-K^y )|v(tu)|^<  2a||v||2,  a = max  (-cr(2nco)c+  R.PC, 
o>=-N  " 0<U|<N  11 


Since  a > 0,  and  if  p > 2,  then  a is  hounded  independent  of  N,  and 
0.22)  and  lemma  3.2  give  us  the  energy  estimate 


(v,v)f  < 2 Real  (v,  (CN  - ^ iiB^/^x)  * v)  + 2a||v||2 


< 2 (max  |CN  - ^ SBN/6x|  + a)V|£L  . 


If  p > 3 then  a simple  calculation  gives  us 


V < < ard  * i)  hi! 


-1  -2 

Therefore,  if  2ttct  > M^K^(N  +N  ) then  o:  in  (3.2?)  is  nonpositive  and  ve 
obtain  the  following  theorem  from  (3*22). 

-1  -2 

Theorem  3*1.  If  P > 3 and  2 TTa  > (M^K^(K  +N  ),  then  the  solutions  of 
(3.1^)  satisfy  the  estimate 


* 


(^.24 ' 


{v>v)t  $ 


2 Feal  lv,  \C. 


- 2 6V6x) 


* v 1 


This  is  entirely  satisfactory  since  it  is  essentially  the  same  as  the 

corresponding  estimate  for  the  differential  equation.  Furthermore,  N can 

1 P 

always  be  chosen  large  enough  so  that  2na>M1K1(N  +N  ),  at  least  in  principle. 

For  hyperbolic  equations,  A = 0,  the  situation  is  not  as  good. 

In  this  case  we  have  to  control  the  smoothness  of  v.  Experience  has 
shown  that  higher  frequency  modes  can  grow  if  this  is  not  done. 

Let  m > 1 be  a natural  number. 


v = 


y ~ frriaac 
v(c«)e 

oj=_n 


and  define  v^,  v0  by 


(3.25) 


vi  = 


Via  ' 


2TTio)X 


\c  | < N1 


v2  = V - V1  ’ 


where  - (1  - l/m)K.  The  smoothing  operator  H = Il(j,m,D)  mapping 

Tjj  into  T.j  is  defined  by 


(3.2b) 


where 


w = Hv  = \ w(ot>)e' 

u>=-N 


27ria)x 


/ " 
v 


(cd)  if  |a)|  < (1  - i)N 


w(oj)  = / v(co)  if  |(o|  > (l-i)N  and  |v(to)|  < 


. D|k||  . 

( 1 


D!|v 


'1 


(2tt|cd|  )<:i 


(27r|ca|  )J  |v(tu)| 


1 1J 


otherwise 


is  a natural  number  and  D is  a constant.  Thus,  only  the  higher 


5 

frequencies  are  modified,  i.e., 


Hvi = vi  > IMI  - IMI  • 


We  want  to  show  that  H is  a very  mild  form  of  smoothing. 

Lemma  3.3.  Let  V > 0 be  a constant  and  j a natural  number, 
the  class  of  functions  with 

(?*27)  Vu/Wl!2  < y2!M!2  . 


(3-28)  (2ttN(^'  )2J  > 2 Y2  and  D y 


then 


Hu  = u 


Proof.  Let  u e and  write  it  in  the  form 


u — u.  + u where  u,  (o>)  = 0 for  |cjd(  > IHli  n 
■L  d 1 m 


V3-27)  implies 


(““■  2nN)2J!!u2!;2  < | |b J" /ft x J' ||2  < y2  ( !!Ul||2  + 


3y  (3-28) 


Consider 
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L 


Therefore,  for  co  / 0, 

|u(co)|2  < Y2(2tt|o)|  )"rj||u||2  < 2y2(2n|co|  )‘2<j||u1||2 
and  the  lemma  follows. 

Instead  of  (3.11)-)  we  now  consider  the  approximation 

(5*29)  3T  = * Hvx  + cn  * v > 

To  see  that  (3.29)  has  a unique  solution  we  need. 

Lemma  3.^.  H is  a Lipschitz  continuous  operator  from  T^  into  T^. 


Proof . Let 

v(i)  0 tn 

and 

w^  = Hv^^,  i - 1,2.  Note  that 

|w(l)(u))i  < 

|v(l)(a))| 

and 

arg  w^(co)  = arg  v^^co),  i = 1,2,  both 

follow  from  the  definition  of  H.  Consider  the  quantities 

ia(1).  . *(  2 )/  t I 

|w  (tuj  - w (od;|.  We  consider  three  cases.  Let 


= {o>|  |to|  < N,  w^\cu)  = v^^(o)),  H - 1,2} 
Jp  = feu | | o)|  < N,  w(i)(o>)  / v^(co),  i = 1,2} 
J3  = [o|  | cd|  < N,  co  f J±  U J2} 


From  the  definition  of  H it  follows  that  u)  e ^ if  |o)|  < = N(  1-1/m). 

If  o)  c J^,  then  |w^\<n)  - w^'  \cd)  | = |v^^(co)  - v^2^(co)|.  If 
o)  t Jn,  then 


d(o))  s 


|w^ 1 ^(03)  — W^  2 ^(o) 


) I = |k( 


CO) 


.(1), 


TTTT^  ‘ K(a)) 


,(2)| 


'Ml 


*(2)/  \ 
7 Mi 
a(  2 ) / 1 

v '(cu)| 


where  K(co)  = D/(2tt  |cd|  )J.  We  assume,  without  loss  of  generality, 

(1)||  ll„(2)l 


that  ||v^  'll  > ||v'1 


Using  the  triangle  inequality  we  obtain 


:(D/ 


\„,  \||  (2)||  V^'(uO  mi  (2)n  V^‘  ^(oo) 

lK(u))Uvi  ll  rd)  , - II  7.(2 r;, 

|vv  '(cd)|  |vv  '(uj)| 


We  can  bound  the  first  term  of  our  last  expression  by 


K(x>)  j ||vp-)||  - Hv^lll  < KCco)!^^1^  - v^||  < K(co)||v^1^  - v^|| 


since  the  two  corqplex  numbers  have  equal  arguments.  We  can  bound  the 

iA(l)/  \ \i 

second  term  by  |v  '(aij  - v'  (co)|  utilizing  the  triangle  inequality 

i®l  i02 

and  the  fact  that  the  distance  between  two  points  re  and  r^e 

J.  d 

is  a non-decreasing  function  of  r^  if  r^  > r^.  Finally,  we  obtain 


(3-30)  d(o>)  < K(oj)||v^1^  - 2 ^ ||  + Iv^^co)  - v^^(cu)| 


if  to  e J0.  Let  o>  e and  assume  without  loss  of  generality  that 


v'  (d)  f-  \r  (d)  and  v''L  (d)  = w^  (d).  If  |v'  '(aa)|  > K(ao)||v^  ||, 

then 

d(o>)  < Iw^^co)  - K(o))||v|1^||  | + |K(a))||v^1^||  ^(2  )|  " \oo)  | 


h(2) 


< |v^(a>)  - \cu)|  + |K(u))||v^1^||  —^2) K(u))||v^'  ^ 


A(  2 ) 

_v__ 

I 2 ) I 
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< Iv^^oo)  - v^(oo)|  + K(o>)  | Hv^1  ^ ||  - ||v[2)||| 

< l 1 ^ ( a) ) - v^(co)|  + K(co)||v(1^  - v^|| 


If  |v^‘  \oo)  l < K(cn)||v|j1)||,  then  it  easily  follows  that 

d(co)  < |v^\a>)  - v^2^(co)|.  Thus,  if  ix>  e J y d(a>)  satisfies  the 

inequality  (3.30).  Now  we  estimate 


« p 

S cl  (to) 
tn=-N 


< 


| v^\cn) 


:(2) 


(co)  |‘ 


coo  J .. 


+ E (K(co)||v(l)-  V(2)||  + 
cdCJ?Uj, 

< (2  + 4K2(N1)(N-N1))||v^1^  - v^ 


i-(l)/  > A(  2 ) , \ | \2 

|vv  ^ (to ) - v (co)|) 


which  yields  the  desired  result. 

From  Lemma  3.4  it  follows  that  the  operator  on  the  right  hand  side 
of  (3.29)  is  Lipschitz  continuous  and  it  then  follows  that  (3.29), 
with  initial  data,  has  a unique  solution  v(t).  We  will  now  derive 
estimates  for  the  norm  of  this  solution. 

We  have 

dJMI2  = 2 Real  (v,vt)  = 2 Real  (v,^  * Hvx  + CN  * v)  . 

The  term  (v,C^  * v)  is  easily  bounded  as  before  using  Lemma  3*2  if 
C f P(o:,m)  with  Ot  > l/2,  or  is  continuously  differentiable. 
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We  write 


* Hvx)  = * (vi)x)  + * ((vi)x  - Hvx) 

splitting  v = v^  + Vg  and  utilizing  the  fact  that  H does  not 
alter  the  first  Fourier  components  of  the  vector  it  operates  on. 

We  then  further  split  ^ * (vp)x  in  terms  of  Q and  R = R^  + Rg 
as  before  to  obtain 


2 Real 


<V'B» 


HVX)  = 


2 Real 


{ ( v>  R1V1 ) 


+ ( v> Rg^ ) + (v,^  * ((v1)x  - Hvx)} 


where  we  have  used  the  fact  the  Q is  skeu-hermitian.  Recall  that 
R-^v^  = - i dB^/dx  * v^  which  is  bounded  as  before  if  B e p(a,M) 
with  a > 3/2.  We  have 


O 1 

(3.31 1 Jt  ' 1 " “ ^ Real  (v,C^  * v - ^ dB^/dx  * v^)  +2  Real  (v^R^v^) 

2 Real  (v,  ^ * ((v^  - Hvx) 


the  first  term  is  bounded  and  converges  to  the  proper  estimate  for  the 
differential  equation.  We  will  now  construct  bounds  for  the  last  two 
terms.  We  assume  that  satisfies  (3.19)  and  obtain,  corresponding 

to  (3*20)  , 


(3.32) 


N N 

| ( v,R0vn  ) | <n(2N+l)(|  r $(«>)  T &..(u)v  (co-?N-l-u)| 

u>=o  m=-h 


-1  N 

+ I z v(u>)  r iI(u)v1(a^2N+l-d)|)  . 

cn=-N  p~-N 


2k 


Utilizing  (3-19)  ve  obtain 


N N 

I e v(cd)  r iUjOv^oj-aj-i-n)!  < 

03=0  |i=-N 


« to-2N-l+N,  o 

E |v(tu)|  E 1|2nMrt3K1(t0-2N-l-^)|  < 

03=1  4=-N 


0-2N-1+N 


Ml  (^j)  £ |v(o>)|  E |v1(u>-ai-i-n)|  < 

co=l  H=-N 


\ (2^)t3N1/2||v1||  E |v(u3)|  < 
03=1 


h <2^  NHVJIM 


and  the  second  term  on  the  right  hand  side  of  (3-32)  also  satisfies  the 
same  estimate.  We  obtain 


13.33) 


l(v,R2v1)|  < 2nM1(2N^+N)(^)S||v1||||v|| 

< (3/(2n)3'1)M1m(3N"P+2||v1||||v| 


We  only  have  the  term  (v,B^  * ((v^)x  ~ Hvx)  left  to  estimate.  We 
have,  via  lemma  3.2,  that 


(3.34)  |(v,^  * ((v1)x  - Hvx)|  < max  | | ||v||||( ( v-]_ )x  - Hvx 


From  the  definition  of  H ve  have 
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1 


l!<vA 


(3.5b) 


D|i(  v ) 

Hv  ) ||  < 2 — 4,  - ■ 

(2-)J  o N +1  | |' 


< 


2D 

(2n)J_1(j-l) 


l|v|| 


if  j > 2. 


We  can  now  collect  our  estimates  (3«5l),  (5.53),  (3*3*0  and  (3.35) 
to  obtain 


Theorem  3.2.  Let  j = 3 > 2,  then  the  solutions  of  (3.29)  satisfy 
the  estimate 


\llv|r  < ? Real  (v,CN  * v - | d^/dx  * v1)  + 


(5.56) 


[(b/(2n);;_1)MLnrr"'  + ( to/(2n  )J"‘L(  J-l)  V"'  max  I^UM! 


If  j - ? > 2,  then  the  estimate  (3.35)  converges  to  the  corresponding 
estimate  for  the  differential  equation  as  N ■*<*>. 

If  the  coefficients  are  smooth  the  estimate  (3*35)  is  quite 
satisfactory  for  sufficiently  large  N.  We  have  been  able  to  obtain 
this  estimate  by  introducing  the  smoothing  operator  H and  by  requiring 
that  the  coefficients  C and  B be  smooth.  A similar  estimate  can 
be  obtained,  with  much  less  effort,  if  we  were  to  alter  the  definition 
of  H such  that  w(o>)  - 0 if  Jco | > N^,  or  w(cn)  - v(co)/(  ( 2n[  |o>|  ]+)J+l ) 

if  |oi|  > where  [g) denotes  the  positive  part  of  g.  These  are 
both  linear  operators.  However,  the  resulting  methods  are  less  accurate. 

Convergence  estimates  can  be  constructed  utilizing  the  estimates 
of  theorems  3.1  and  3.2  following  those  of  Kreiss  and  Oliger  If]  and 
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Fornberg  [5J  and  the  approximation  results  of  Bube  [1], 
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